Effect of electron-phonon interaction on spectroscopies in graphene 
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We calculate the effect of the electron-phonon interaction on the electronic density of states 
(DOS), the quasiparticle properties and on the optical conductivity of graphene. In metals with 
DOS constant on the scale of phonon energies, the electron-phonon renormalizations drop out of 
the dressed DOS, however, due to the Dirac nature of the electron dynamics in graphene, the band 
DOS is linear in energy and phonon structures remain, which can be emphasized by taking an 
energy derivative. There is a shift in the chemical potential and in the position in energy of the 
Dirac point. Also, the DOS can be changed from a linear dependence out of value zero at the Dirac 
point to quadratic out of a finite value. The optical scattering rate 1/r sets the energy scale for 
the rise of the optical conductivity from its universal DC value Ae'^/irh (expected in the simplest 
theory when chemical potential and temperature are both <C 1/2t) to its universal AC background 
value ((To = -Ke^ /2h). As in ordinary metals the DC conductivity remains unrenormalized while its 
AC value is changed. The optical spectral weight under the intraband Drude is reduced by a mass 
renormalization factor as is the effective scattering rate. Optical weight is transferred to an Holstein 
phonon-assisted side band. Due to Pauli blocking the interband transitions are sharply suppressed, 
but also nearly constant, below twice the value of renormalized chemical potential and also exhibit 
a phonon-assisted contribution. The universal background conductivity is reduced below ctq at large 
energies. 

PACS numbers: 78.67.-n,71.38.Cn,73.40.Gk,81.05.Uw 



I. INTRODUCTION 

Graphene consists of a monolayer of a honeycomb lat- 
tice of carbon atoms. It has been studied theoretically 
since the early work of Wallaceii An important aspect of 
the charge dynamics in graphene is that it is governed by 
a Dirac rather than a Schrodinger equation. Fcrmionic 
energies are proportional to momentum with the effec- 
tive speed of light vq ~ lO^m/s. While the material was 
isolated experimentally only in 2004;^ many fascinating 
discoveries have already appeared^i^ including an inte- 
ger quantum Hall effect with half integer filling factors, a 
minimum conductivity and a Berry phase shift of TTi^i^!^ 

An important dimension in the study of graphene is 
that it can be incorporated into a field effect device^ and 
doped away from half filling by changing the gate volt- 
age. This means that the position of the chemical poten- 
tial can be varied with respect to the Dirac point where 
conical valence and conduction band meet at two specific 
K points in the Brillouin zone. As is the case in ordi- 
nary metals the electron-phonon interaction^"'^^ (EPI) 
renormalizes the bare bands and changes the properties 
of graphene. Here we consider specifically the electronic 
density of states (DOS). In metals, where the DOS is 
essentially energy independent on the scale of phonon 
energies, a well-known result is that the electron-phonon 
self-energy entirely drops out of the problem and the DOS 



retains its bare value. In graphene, as we will show in this 
paper, this is no longer the case and there is an imprint 
of the EPI in the renormalized DOS. This is expected in 
systems where the DOS varies on an energy scale com- 
parable to phonon energies. It also arises in finite band 
systems, with the top and the bottom of the band are 
modified in a particularly important way by interactions 
with the phononsJ^ii^ii^ii^ 

Another known result in simple metals is that the EPI 
renormalizationsiSiiiiii drop out of the DC value of the 
conductivity but significantly renormalize its AC value 
and phonon-assisted Holstein processes become possible. 
For graphene the DC conductivity which is due to both 
intra and interband transitions is also unaffected by the 
EPI. However, the intraband AC Drude contribution has 
its optical spectral weight reduced by a factor of (1-|-A°*^), 
where A''''^ is the effective mass renormalization parame- 
ter (dependent on value of the chemical potential), and 
its scattering rate is also reduced by the same factor. 
Optical spectral weight is shifted to phonon-assisted Hol- 
stein side bands originating from the incoherent part of 
the electron spectral density. The interband transitions, 
which are strongly suppressed by Pauli blocking below 
twice the renormalized chemical potential, also exhibit 
phonon-assisted side bands which add on to the intra- 
band Holstein processes at photon energies above the 
phonon energy. The universal value of the background 
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conductivity at higher energies is suppressed below its 
bare band value ctq — 7re^/2/i. 

In Section II, we introduce the model of Park et 
g^^ ^is, 19^20 fpj. ^j^g electron-phonon interaction in graphene 
which we write in a more general form so as to include the 
self-consistent electronic density of states (DOS) in the 
self-energy itself. Iterations of the self-energy to obtain 
a self-consistent value of the DOS is known to be im- 
portant in finite bands particularly in the region of the 
band edge. The electron-phonon interaction renormal- 
izes the electronic motion at the Fermi surface through a 
renormalization factor (1 -f A''*).^-'^ In addition, the chem- 
ical potential is changed from its bare band value as is 
the relative position of the Dirac point. In the limit of 
small chemical potential, analytic formulas are obtained 
for these renormalizations and are compared with more 
general numerical results. Our results for the renormal- 
ized DOS show phonon signatures at small energies and 
a renormalization of the band edge with renormalized 
bands no longer ending abruptly as in the bare band. 
An analytic formula is derived for the modification of 
the Dirac point due to interactions. 

In Section III, we consider more specifically phonon 
structure in the renormalized DOS, a topic on which 
we have some preliminary and complementary work.^^ 
To this end, we introduce a Lorentzian model for the 
phonon distribution instead of the simple Einstein model 
of Park et al.^^ This provides a more realistic descrip- 
tion of the structure that is expected to be seen in 
experiment 1^ To see this structure more clearly, we take 
a first derivative and find that the Lorentzian line shape 
of our assumed boson spectrum is faithfully reproduced 
in these plots except for a small shift in energy. Sec- 
tion IV deals with aspects of quasiparticle renormal- 
ization as they show up in angle resolved photoemis- 
sion spectra (ARPES).— i^^i^^i^l We particularly empha- 
size the shift in chemical potential^^ revealed in such 
curves and note that electron-phonon renormalizations 
remain to high energies. Section V deals with optical 
conductivity2ai^i^i^i2^iMi2^i^i21i2a starting with the ef- 
fect of the electron-phonon renormalization on the low 
frequency conductivity. We begin with the DC conduc- 
tivity which is found to be unrenormalized. Conditions 
for the observation of the universal limit are discussed 
and the contribution to this value of intra and interband 
transitions are separately determined. The rise of the 
conductivity cr(f2) from its universal DC to its universal 
AC background valued with increasing frequency is de- 
scribed. Analytic formulas are provided for the electron- 
phonon renormalization that arises for finite doping when 
fi is much larger than the impurity scattering rate 1/r. 
We find that in this case, the intraband transitions pro- 
vide a Drude-like contribution to the conductivity which, 
just as is known to be the case in ordinary metalsjii is 
renormalized by the mass enhancement factor (l-t-A***^) in 
two ways. The effective plasma frequency of the Drude as 
compared to its free band case is reduced by 1/(1 + 
and the scattering rate l/r also goes into {1/t)/{1 + X'^^). 



In addition, optical spectral weight is transferred into 
Holstein phonon- assisted absorption sidebands. On the 
other hand and in sharp contrast, the interband transi- 
tions at small Q remain unrenormalized but beyond the 
phonon energy also exhibit Holstein sidebands. In sec- 
tion VII, the infrared region of the conductivity^i^i^i^ 
is considered more explicitly with the specific aim of un- 
derstanding separately the role played by intra and inter- 
band transitions. It is explained how interband transi- 
tions which would be Pauli blocked in the pure bare band 
case can nevertheless take place in the interacting system 
although with reduced spectral weight. This arises be- 
cause in the interacting system Bloch states are never 
occupied with probability one and so Pauli blocking is 
not complete. The partial sum rule on the optical con- 
ductivity is discussed. Section VIII contains a summary 
and conclusions. 



II. THEORETICAL BACKGROUND 

In graphene, the massless Dirac nature of the electrons 
means that their energy (et) is linear in momentum (Kk) 
and consists of two branches, i.e., ek = ±?i?;o|k| with the 
velocity vq equal to about 10^ m/sec and k the wave 
vector measured from one of the two K points in the 
Brillouin zone where valence and conduction bands meet. 
The ± corresponds to the upper and lower Dirac cones or 
particles and holes, respectively. Here, e can serve as the 
label for the absolute value of momentum. The carrier 
spectral density A(e, uj) can be written in terms of the 
self-energy S(w) as 
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-Iml](a;) 



TT - Rel](w) + fi-e]^ + [Iml](w)]2 ' 



(1) 



Valence and conduction cones are both part of Eq. ([J) 
and correspond, respectively, to negative and positive e 
label. As defined, S(cj) is the self-energy such that lu = 
corresponds to the Fermi level and fi is the chemical po- 
tential of the interacting system. The electron-phonon 
interaction in graphene has been discussed by many 
authorsi^ii^i^^ including recent work^2i^i2i directed to- 
wards understanding its effect on the optical conductiv- 
ity. Here we find it convenient to follow the formulation 
of C.H. Park et alM- later also generalized to bilayers 
and graphitCi^-^ These authors proceed within an ab initio 
pseudopotential density functional calculation of the elec- 
tronic bands in the local density approximation. They 
further calculate the phonon frequencies and polariza- 
tion as well as the electron-phonon matrix elements from 
which they construct the electron self-energy. The impor- 
tance of this work to the present discussion is that in the 
end they provide a simplified model for the self-energy 
which they find captures all essential elements of their 
sophisticated numerical calculations. They show that to 
a good approximation, one can think of the electrons 
as coupled to a single phonon of frequency uje = 200 
meV and provide a value for the relevant coupling. The 
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self-energy is independent of electron momentum label 
as well as band index, electrons or holes, a simplification 
they tracei^ to the Dirac nature of the electronic states 



in graphene. Based on these simplifying ideas, we begin 
with a self-energy for temperature T = of the forniHiii 



, ,N(uj') A 



No Wc 



e{j) 



u) — uj' — uje + «0+ uj — uj' + uje + iO~^ 



(2) 



In Eq. ([2]), lue is the Einstein oscillator energy, uj is in 
units of energy, A is the coupling, 0{lo) is the Heaviside 
function, and No = 2/i:%^Vq. Wc is the cut off on the 

Dirac cone given by \/ 7r-\/3i, with t the nearest neighbor 
hopping parameter taken to be 3 eV which corresponds 
to Wc ^ 7000 meV. In Eq. we have added an addi- 
tional element not included in the previous work of C.H. 
Park et alj^ but which can be important in the case of 



graphene: N(uj) is the self- consistent electronic density 
of states given by 



Nju;) _ \e\ ^ImI](co) 

No ~ J.wc ^ ^ - ReS(c^) + /i - e]2 + [ImE(w)]2 ' 

(3) 

This can be evaluated to give 



N{uj) 
No 



arctani - 



Wc 



' arctan 



Wc 



— 2 arctan 



27r 



■In 



Wcf + T^W 



(j)2 + r2)2 



ill 



where T = — ImE(aj) and Co — lj — Rel](aj) + ^. In 
good metals, the electronic density of states (DOS) in the 
energy range important for the electron-phonon interac- 
tion around the Fermi energy, does not vary significantly. 
Only its value at the chemical potential is relevant. In 
such systems, the linear in e factor |e| is not present in 
Eq. ([3]) and extending the integration limits to infinity 
gives a constant independent of ImS(a;) and ReS(u;). In 
more complicated metals such as the A15 compounds, 
as an example, it was recognized early on that such an 
approximation is no longer validi2, and in this instance it 
is the self-consistent DOSHiiiii^ that enters Eq. ©, as 
in the work of Engelsberg and Schrieffer;^^ This is also 
relevant for finite bands i^'^i^^'^^ For graphene, the bare 
band structure density of states is linear in energy rather 
than constant and hence we expect significant changes 
could be introduced into when equations ([2]) and 

([3]) are iterated to convergence. We will test this in what 
follows. 

As a first approximation, we can replace N{lu) in 
Eq. ([2]) by its bare non-interacting value given by 

N{uj) _ f |w + /zol, for - VFc - Mo < w < VFc - Mo, 
A^o 1 0, otherwise, 

(5) 



where mo is the chemical potential for the bare bands, 
i.e. without electron-phonon rcnormalization. All states 
with UJ < are occupied and with a; > unoccu- 
pied, so that UJ — sets the boundary between occu- 
pied and unoccupied energies and by choice this remains 
the case even when interactions are included. Note that 
the value of mo sets the doping p (per unit area) with 
p = sgn{fio) fJ-o / {TTJi^ Vq) , with mo negative for holes and 
positive for electrons. This is illustrated in the bottom 
frame of Fig. [1] by the black dotted lines for the case of 
ImoI = 150 meV. The Dirac point is at -150 meV (-^150 
meV) for electron (hole) doping and corresponds to zero 
DOS 01 UJ = —po from Eq. ([5]). For this bare N{uj) the 
self-energy of Eq. ^ can be evaluated analytically but 
is different for po > and po < 0. For po > (writing 
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ReS 



Mo>0 



(CO) 



^/ 1 



{Wc + UJE - - ^J'o)i^J■o +UJ + loeY 



-(/xo + w)ln 



- ul.){Wc + + t^_E + Mo) 



(6) 



r 



and for negative /io 

ReE^o<o(w) 



Rel]|^„l(-w). 



(7) 



This symmetry between positive and negative values of 
chemical potential also holds in the case when the full 
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self-consistent density of states (DOS) N{u)') is used in 
Eq. ([2]). Before showing this explicitly, we note that 



Iml](w) 



■——\LO-LOE+^iQ\, for LJE < <Wc - + OJE, 

Wn 

IT J± 

— — |cj + cjb + /io|, for - ljb > ^ > -VFc - Mo - ^^B, 
I Wc 



(8) 



r 



and zero outside these intervals. For negative values of 
Mo we note 



ImS^o<o(w) = ImE|^„|(-w). 



(9) 



Using Eqs. ([7]) and (O in Eq. ^ for the self-consistent 
DOS, we see that 



iV^„<o(c^) = A^|p„|(-t^) 



(10) 



and this relationship implies that the symmetries of 
Eqs. ([7|) and ([9|) apply even when the full self-consistent 
DOS is used in Eq. ^ for the self-energy. 

In Fig. [1] we show results for ReE(ijj) (top frame), 
Iml]((jj) (middle frame) and N{lo)/No (bottom frame), 
where Niiu) is a first iteration of Eq. ^ which should 
be sufficient for small enough electron-phonon interac- 
tion. The dashed (blue) lines are for mo > (electron 
doping) and the solid (red) for mo < (hole doping). 
The dotted (black) line in the bottom frame is the bare 
band density of states shown for comparison. Here the 
chemical potential mo was taken to be |mo| = 150 meV, 
Wc = 7000 meV and cj^ was taken to be 200 meV as 
in Ref.^*. The value of the electron-phonon A parame- 
ter was set at 250 meV. The corresponding spectral mass 
enhancement A = 2^4/0^^; ~ 2.5. As we will show, this pa- 
rameter is very different from the true mass enhancement 
which we will call A''^ and which describes for graphene 
the renormalization of the velocity because the carriers 
are massless. 

Note that the symmetry incorporated into Eq. (O 
is satisfied in the top frame for ReE(a;). We see 
logarithmic-type singularities^'^ at w = ±W£; but these 



would be smeared if one had used a distributed phonon 
spectra rather than an Einstein model. The logarithmic 
singularity can be traced to a factor {uj ±loe) in the ar- 
gument of the logarithm in Eq. ^ and these will lead 
to sharp structure in the electronic density of states. For 
later reference, we note an important exception. For zero 
chemical potential mo = 0; cancellations occur in Eq. ^ 
and the singularities are removed. 

ReE(a; = 0) is not zero but -24.7 (-^24.7) meV for 
Mo > (mo < 0). This quantity relates interacting (m) 
and noninteracting (mo) chemical potentia l^^i^^i^^ with 



M = Mo + ReE(w = 0). 



(11) 



If for simplicity we assume all three energy scales w, Mj 
and loe to be small as compared with the large band cut 
off Wc, then Eq. ^ for ReS(w) gives at w = (mo > 0) 



2A ( 

ReE(cj = 0) = — <^cjBln 

Wc L 



Mo + 



UJE 



- Mo In 



Wc 



Mo + (^E 



(12) 

The shift between interacting and noninteracting chemi- 
cal potential is zero only for the case mo = (particle-hole 
symmetry). For <ti uje 



ReY.{uj = 0) 



2A 



Mo 



In 



1 



(13) 



In Eq. p^ . Wc is involved, the limit Wc oo can- 
not be taken, and ReS(Lj = 0) is a fraction of mo be- 
cause A/Wc is a small number. In terms of the spec- 
tral A introduced earher (A = 2A/(jJe), the coefficient 
2A/Wc = X{uje/Wc) with lue/Wc ~ 1/35 for realistic 
parameters corresponding to graphene. 
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FIG. 1: (Color online) Real (top) and imaginary (middle) part 
of the self-energy E(w) and the density of states N{u!)/No 
(bottom) (all in units of meV) as a function of u} in eV for 
1^0 1 =150 meV. Dashed is for fio > and solid is for /io < 0. 
The bare band DOS is indicated by the dotted curve. 



For small but finite values of lo, we get from Eq. 
in the same limits, (|/zo| < uje), 



Rel](w) 



2A 



In 



Wc 



Mo + LOE 



Imp I 



l-HZli I +ReI](w = 0). 

(14) 

Denoting the carrier effective mass renormalization due 
to the electron-phonon interaction by then for uj 0, 
the self-energy is given by 



Rel](w) = -X'^'^Lu + ReS(w = 0). 



(15) 



The renormalized energy (Ek) as measured, for example, 
in angle-resolved photoemission spectroscopy (ARPES) 
is related to the bare energy near the Fermi energy by 



Ek = 



1 + X 



off 



(16) 



with hkF, the Fermi momentum of the bare band, so 
/io = -izhv^kp, with the -I- sign for electron doping and 



— for holes. We see that interactions renormalize the 
bare effective speed of light wq to Vq = vq/{1 + \°^). 
To see this, we return to the spectral density of Eq. ([1]) 
and note that bare and renormalized energies near the 
Fermi surface are obtained when the first bracket in the 
denominator is set equal to zero: 



Ek + X^^Ek - RcJ:{uj = 0) + m = fife, 



which gives Eq. 
we get 



(17) 



16]) with (fc — kp) small. For /ig ^ i-^e 



.off _ 2A 



In 



c 



LUE 



(18) 



which is ^ 0.19 for the parameters of Fig.[T] This value 
is larger than calculated in density functional theory — 
but smaller than measured experimentally^''^. Note that 
in terms of Eq. (US]), Eq. (fT^]) can be written as 



RcS(tj = 0) = -MoA' 



cff 



and so noting Eq. pT|) 

M = Mo(l - A^*^) 



Mo 



1 + A«ff 



(19) 



(20) 



since A'^*^ is small so that in this simple limit the factor 
(1 -I- A°*) also renormalizes the value of the bare chemical 
potential. 

The middle frame of Fig. [1] gives first iteration results 
for the imaginary part (ImE(a')) of the self-energy which 
clearly satisfy the symmetry of Eq. Iml](w) is zero 
between —uje < lu < uje at which point it jumps up to a 
finite value because only then can a quasiparticle decay 
by boson emission. In an ordinary metal the scattering 
would remain constant above lue because the final den- 
sity of electronic states aX lj — uje is constant. This is not 
true for graphene where N{lo~loe) varies and the scatter- 
ing rate reflects the energy dependence of the underlying 
bare band structure seen in the lower frame of Fig. [1] as 
the black dotted curve. Looking at the red solid curve 
for energy less than wg, the |ImE(aj)| increases linearly 
with uj as does iVo(a;) and for w > we see first a drop 
towards the zero in |ImI](a;)| after which it too increases 
linearly, imaging the bare density of states as captured 
in Eq. ([5]). In more reahstic calculations there would be 
a nonzero contribution to the self-energy from acoustic 
phonons and 1tsiY,[uj) would never be exactly zero as we 
have it here. But from the work of Ref J^, this is small in 
graphene. Finally, we note that for a fixed value of ±/io, 
ImE(w) is not symmetric with respect to w ^ —oj as it 
would be in good metals. Here such a symmetry results 
only when /io = 0, i.e. the case of no doping. 

In the lower frame of Fig. [TJ we show our one iteration 
results for the DOS A^(a;) (solid red curve for /io < and 
dashed blue for /io > 0) and compare with the bare case 
(dotted black curves). N{uj) is not symmetric about lu = 
and the renormalizations are not the same for electron 
and hole branches. Also, the Dirac point is shifted to 
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higher energy as compared to its bare band position for 
/Xq > and to lower energy for /xg < 0. The shift becomes 
greater with increased doping. Another feature to notice 
is that bare and renormahzcd DOS have exactly the same 
value at a; = 0, which is where by choice we have taken 
the Fermi level to be in both cases. This makes sense 
because the imaginary part of T.(lu) at o; = is always 
zero and so the Lorentzian forms in Eq. ([3]) reduce to 
delta functions and for fiQ > 0, as an example, 



N{uj = 0) 



^ — Rel](a; = 0) = ^o, 



(21) 



which is the value of the bare density of states at the 
bare chemical potential. The value of the DOS at the 
Fermi surface remains pinned to its noninteracting value. 
Another important feature of our DOS results is that 
phonon structures are clearly seen as small kinks in the 
curve. This will be smeared somewhat when a distributed 
spectrum is used in computing the self-energy rather than 
the Einstein mode of Eq. We will return to this is- 
sue in a later section where we focus more specifically on 
phonon structures. For the moment, we note that our 
results for N{cij) versus oj are in striking contrast to what 
is found to apply in conventional metals. In that case, 
the density of electronics states around the Fermi surface 
is constant on the scale of a few times the phonon energy 
and no renormalization of the DOS due to the electron- 
phonon interaction is expected or observed. Further in 
conventional metals, the energy bands are rigidly filled as 
the chemical potential is varied. For graphene, when cou- 
pling to a boson spectrum is included, the bands become 
renormalized differently for each value of the chemical 
potential /ip. These are not rigid bands. 

Next we look at how the position of the Dirac point 
shifts from its bare value at w = —fj-o when interactions 
are included. In the interacting system, the Dirac point 
can still be identified with zero momentum fc = (mo- 
mentum label) which corresponds to e = (energy la- 
bel) as e = ±?iuo|k|. If there was no broadening (but 
of course there is), the spectral function would become 
infinite (peaks) at w — ReE(cj) + /_* = 0. If we write 



ReE(w) = Rel](w = 0) + ReST.{Lo) 



(22) 



where Re5I](a') = at a; = by arrangement, we get 

u) - ReE(w = 0) + A* - RcSJ:{u;) = (23) 

as the transcendental equation which needs to be solved 
to obtain the energy of the Dirac point in the interacting 
band. For one iteration, ReS(ti;) is given by Eq. ^ from 
which Re(5I](tj) can easily be extracted. Noting Eq. pT|) . 
Eq. (QSl) simphfies to 



(24) 



where LOd is the energy of the Dirac point in the interact- 
ing system. For the bare band case KedJ^{uJd) in Eq- (EH) 
is zero and we recover the known result that the Dirac 
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FIG. 2; (Color online) Results of a first iteration for the shift 
in chemical potential ReE(a) = 0) (top frame), the shift in 
position of the Dirac point Re5E(wti) (middle frame) and the 
imaginary part of the self-energy lmS{uJd) at the Dirac point 
(lower frame). All of these quantities are shown in units of 
meV. 



point is at — /io- Re(5S(cL'rf) which we plot as a function 
of /io in the middle frame of Fig. [2] gives the shift in ujd 
away from — /zq resulting from the electron-phonon inter- 
action. It has the same sign as fio and hence displaces 
the Dirac point to the right of its bare value for fj,o > 
and to the left for /io < as we noted before. In the 
upper frame of Fig. [2] we show a related quantity namely 
Rel](a; = 0) which is the shift in chemical potential (/i) 
from its bare band value (/io)- It carried the opposite 
sign to /io and so shifts /i to the left of /xq for positive 
fiQ and vice versa for negative fiQ. This shift in chemi- 
cal potential is expected because the interactions modify 
the density of states and hence to keep the carrier im- 
balance fixed, the position of the Fermi level needs to 
be modified. The question of how this shift might man- 
ifest itself will be addressed in a later section. Shifts in 
the Dirac point position can be seen in ARPES^^ and in 
scanning tunneling microscopy (STM)'*'^. Noting that the 
electron-phonon interaction makes only a small contribu- 
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tion to the measured total quasiparticle scattering rates2£ 
seen in graphene and that the electron-electron interac- 
tions are needed to understand these, i.e., particle-hole 
and plasmons,^^ any conclusive comparison with experi- 
ment would need to include these additional interactions 
and so we do not attempt here a comparison with data. 
Before leaving this section, we note that for /iq ^ ^^e we 
can show that 



200 



Mo 



l-l-A'' 



(25) 



from ([24| with KeSY^{LLjd) replaced by Re(5I](— /io)- Again, 
to a first approximation, the (1-1- A''*^) factor renormal- 
izes the position in energy of the Dirac point. In this 
approximation, the magnitude of lu^, i.e., \LOd\, is equal 
to the magnitude of the renormalized chemical potential 
Except for prominent phonon structures seen in the 
middle of Fig[2l not present in the top frame, |Re(5S(a;d)| 
and |ReI](a; = 0)| are not very different in magnitude as 
compared with the magnitude of /i itself and so renormal- 
ized values of \ujd\ and are close to each other, even 
for a general value of chemical potential. 

In Fig. [3l we consider the case of /io = 500 meV and 
compare our initial step uniterated results shown as the 
dashed lines with results obtained (solid curves) when 
DOS (Eq. ^) and self-energy (Eq. Q) are iterated to 
convergence. We see some changes in real and imaginary 
part of 2(0;) top and middle frames, respectively. We 
note however that the change in slope of the ReS(ti;) out 
of zero is not changed much in the solid curve which 
means that the one iteration estimate of the effective 
mass renormalization A°*^ is quite good. Results for X"^ 
are shown in Fig.[?]as a function of ^q. Its value increases 
by a factor of two over the range shown. 

The differences between iterated and noniterated re- 
sults for the —lmY,{uj) shown in the middle frame of Fig. [3] 
are larger and can become significant. It is clear that 
self-consistency could be important in some quantitative 
estimation of electron-phonon effects on the quasiparticle 
scattering rate and that LDA can underestimate these in 
some cases. There are also changes in the resulting DOS 
shown in the bottom frame but these are small in the 
case shown. Effects of self-consistency become essential, 
however, when the band edge is considered. Fig. [5] is the 
same as Fig. [3] but now the entire band width is shown 
within the linearized Dirac cone approximation with cut 
off Wc- On this scale the phonon structures visible as 
sharp peaks in Fig. [3] are hardly seen. Note that in the 
one iteration case the imaginary part of S(ci;) (middle 
frame) dashed curve becomes zero at the bare band edge 
because the DOS runs out at this point but the solid 
curve which results when we iterate remains finite outside 
this range before eventually going to zero. These features 
of the renormalization get reflected in the density of state 
shown in the lower frame of Fig. [5] in which we have also 
added for reference the bare DOS No{uj) as the dotted 
(black) curve. It is clear that the electron-phonon in- 
teraction profoundly renormalizes the band edge. These 
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FIG. 3: (Color online) Same quantities as for Fig. [T] but 
now the self-energy and density of states have been self- 
consistently iterated and /io = 500 meV which is greater than 
uje- The solid curves are for the iterated case and the dashed 
curves are for the initial uniterated results. The dotted curve 
is the bare density of states. All quantities on the y-axis are 
in meV. 



effects have been studied before for other models of fi- 
nite bands in, for example, RefSfi^ andi^. To obtain a 
good characterization of the renormalized band profile 
(red curve) around the bare band edge Wc, it is essen- 
tial to iteratcji^iii^ We note that the top of the renormal- 
ized band extends to higher energies as compared with 
the bare band and the bottom extends to lower ener- 
gies. It is the phonon energy which sets the scale for 
this smearing beyond the bare band edge. Of course the 
total number of states must be preserved, i.e. the area 
under N{llj)/No is the same for bare (dotted black) and 
renormalized (solid red) curves. Due to interactions, i.e. 
damping, the states below Wc get depleted and so must 
reappear at higher energies leading to an increase in the 
effective band width. We have verified that to within our 
numerical accuracy the sum rule on N{uj)/No was indeed 
satisfied. As the bands distort, and these distortions are 
different for every choice of chemical potential ^0, it is 
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FIG. 4: (Color online) The mass renormalization parameter 
A"*^ at the Fermi level as a function of bare chemical potential 
Mo- 



clear that the Fermi level will also need to be shifted. 
Again, to within the accuracy of our numerical work, we 
have verified in a few instances that integration of the 
occupied part of the states up to ^jl given by Eq. pT|) 
does give the right charge imbalance p. 

Returning to uj of order, at most, a few times loe, a 
different point is illustrated in Fig. [51 where we compare 
results for the DOS for the case of Fig. [T]with /iq < oje 
(top frame) with a case in which all parameters are kept 
the same but a value of /iq > namely /ig — 500meV 
is used. In the top frame the DOS at the Dirac point 
remains zero because for /zq = 150meV, the imaginary 
part of is still zero at the energy of the Dirac point 
as is shown in the bottom frame of Fig. [5J But in the 
bottom frame for /ip = 500 meV, this is no longer the 
case and the Dirac point becomes smeared and N[ijjd) is 
finite. Note that the shift between Dirac point position 
for bare and interacting bands has increased as compared 
to the case of the top frame (see middle frame of Fig. [2]). 
It is useful to get a simple analytic expression for the 
DOS around the Dirac point. 

For UJ near uod and e near e = 0, the spectral functions 
that determine the DOS are 



J 



-Iml](cj<j) 



ReS(wd) 



ReS]'(wrf)(w 
I 



(26) 



Denoting 1 — JicYi'iuJd) by Z, where S'(wd) 
^\uj—ujd7 cind — IniZj(iU(/) by P, we obtain 

1 r 



and approximately (for \{lu — LUd)Z\ <C ^)r^ 



Njco) _ 2r 



In 



Wc 



r 



ttF 



(27) 



(28) 



where the first term is the value right at uj — ujd and the 
UJ variation off the Dirac point is quadratic in lj — ujd- 
Note that N{u;d) ^ as F ^ in Eq. Also from 

Eq. (HZI), as r ^ 0, A{±e,uj) - 6[{lj ~ ujd)Z ± e] which 
gives for /iq < a linear variation of N(ijj) out of ujd 
with slope modified by Z . We have tested these analytic 
results against numerical work and offer a comparison in 
the bottom frame of Fig. [SI The dashed red curve fits 
very well the solid blue curve in the region of the Dirac 
point. To end this section, wc comment on the range 
of validity of the approximate but very useful Eq. ([^0]) 
and ij^S]) for renormalized chemical potential and Dirac 
point position, respectively. In Fig. [3 we compare exact 
results for p, (solid red curve) and ujd (long-dashed blue 
curve) with the approximation ([^0)1 and , respectively 
(i.e., Aio/(l + A'=*), short-black curve). At small /iq, the 



agreement is excellent as we expected. The deviations at 
higher values of ^ 350 meV should be noted and could 
be of importance in some applications. Finally note that 
on the scale shown, ji and \u!d\ track each other well. 



III. MORE REALISTIC PHONON SPECTRA 

In a metal with constant density of states in the im- 
portant energy range for phonons, the electron-phonon 
interaction does not renormalize its value.-- The essen- 
tial argument can be seen from Eq. ^ which is to be 
integrated over energy, but N{{))deA{e,Ljj) ~ N{0) 
independent of uj. This means that phonons will not 
show up in normal density of states spectroscopy. But in 
graphene, this no longer holds for two reasons. First, the 
Dirac density of states is linear in e rather than constant 
as can be seen in Eq. ^ and after integration S remains 
in formula ([4]). Second, we have finite bands and it has 
been noticed by Knigavko et al>^ that this also provides 
a mechanism whereby an image of the phonon structure 
appears in the electronic density of states. To examine 
phonon structure it is convenient to move away from cou- 
pling to a single Einstein phonon as we have done so far in 
Eq. ([2) for the self-energy. For a distribution of phonon 
energies P{i^), it is necessary to average Eq. ^ over the 
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FIG. 5: (Color online) Same as for Fig. [3] but shown over a 
larger range of lj spanning the band edge. 



desired distribution. As a model, we use the truncated 
Lorentzian form of Refii^: 



P{y) 



A' 



e{5,-\LOo-y\), 
(29) 

where it is peaked at loq with width 5 and truncated at 
energy ± (5c. A' is adjusted to give A used previously. 
With this distribution 



Eior(t^) = / P{v)Y.{uj,v)dv, 



(30) 



where Yi{ijj,v) = Yi{lo) of Eq. ^ and loe of that equa- 
tion now becomes a variable v. Also, the A of Eq. ([2]) is 
omitted in favor of Eq. (P^)) which enters into Eq. ([50)1 . 
Sior(w) is then used for the self-energy in the DOS cal- 
culation using Eq. In our calculations, we have used 
LOo = 200 meV, 5 = lb meV, 5c = 30 meV. 

To see better the phonon structures encoded in the 
electronic density of states N{ijj) of Figs. [T] and[3l it is 
convenient to take a first derivative. This is shown in 
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FIG. 6: (Color online) The density of states N{uo) (solid line) 
vs. Lo for UE = 200 meV, A = 2.5, Wc = 7000 meV. The 
top frame is for = 150 meV < loe and the bottom for 
/io = 500 meV > u)e- The dotted curve is the bare band case. 
For /io > ti^B, N{uj) at the Dirac point is nonzero and becomes 
quadratic. A comparison of this approximate quadratic be- 
havior given in Eq. ((28} is shown as the dashed (red) curve. 



Fig. [5] for several cases. d[N{uj)/No]/du! is normalized 
such that in the bare band, it is one at energies above 
the Dirac point and minus one below, with a jump at 
that point. This is represented in Fig. [8] as the black 
dashed curves with fiQ — 150 meV for the top and middle 
frames and /io = 500 meV in the bottom frame. The 
solid blue line contains the phonons. In the top frame, 
we see a prominent phonon structure at a; = iij^o in 
d[N{uj)/No]/du} which is superimposed on the bare band 
background above which it rises by more than a factor of 
2. On the negative energy side there is a further smaller 
structure at (— ^o~'^o)- AH three phonon structures can 
be traced to logarithmic singularities in the self-energy 
of Eq. ^ of the form \n\u! ± uJo\ for the first two and 
{u! + ^0 + loq) In \u! + ijLq + loqI for the last which is a 
weaker singularity. Note that the vertical drop in the 
solid curve is not at lu = — /ig but is rather shifted to the 
right because the Dirac point has been shifted and this 
drop signals the position of ujd rather than of —fj-o- The 
value of this derivative a,t u; — has height equal to (1 -t- 
A''*^) and this can be used to measure this renormalization 
parameter. Of course this is only possible if one knows 
the value of the bare band Fermi velocity vq. 

On the scale of this figure, it is clear that the phonon 
structure should be easily detectable. In the inset indi- 
cated by the solid arrow, we compare the structure in the 
first derivative of the DOS (blue solid curve) with the in- 
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FIG. 7: (Color online) Comparison of ^ (red solid curve) 
and \ijjd\ (long-dashed blue curve) with approximate relation 
Mo/(l + -^°^) (short-dashed black curve) as a function of /lo- 
The dotted line is for comparison and gives unrenormal- 
ized. 



put spectral density (red dashed curve) scaled down by 
a factor of 10. Except for a small shift between the two 
curves, there is excellent agreement and conclude that 
the first derivative not only indicates where the boson 
structure lies but also captures well the correct profile of 
the Lorentzian spectra that we have used. This is further 
emphasized in the middle frame where two Lorentzians at 
sHghtly different energies wqi = 150 meV and wo2 = 200 
meV are used with the second peak height chosen to 
be twice that of the first. Again the reproduction of 
the details of the input electron-phonon spectral density, 
which is scaled down by a factor of 10, are excellent. 
This demonstrates that normal density of states spec- 
troscopy can potentially be used to probe the phonon 
structures in graphene. The technique could also be ex- 
tended to other structures such as those due to hole- 
particle excitations and plasmons. The lower frame is 
the same as the top frame but now /ig = 500 meV be- 
yond the phonon energies of ~ 200 meV in our model. 
In this case the phonon structures are even bigger as 
we expect since the electron-phonon interaction reflects 
the bare density of states which increases with increasing 
energy. This expected increase in coupling can clearly 
be probed in density of states measurements. A feature 
to be noted is that because /j,o > wq in this frame, the 
Dirac point is changed to the quadratic behavior shown 
in the lower frame of Fig. [H] and this leads to a drop 
in dN(uj)/duj which is no longer vertical but instead is 
smeared. Finally, we note that while the first derivative 
already gives a rather good image of the shape of the 
underlying electron-phonon spectral function, a formal 
inversion of Eq. ^ along the lines used in superconduc- 
tivity by McMillan and Rowell^ could be used to get its 
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FIG. 8: (Color online) Derivative of the electronic density of 
states. Dashed black line is the bare case and the solid blue 
curve includes phonons. The top frame is for a bare chemical 
potential fio = 150 meV and coupling to a Lorentzian phonon 
distribution peaked at uo = 200 meV, while the middle is the 
same but a distribution with two Lorentzian peaks and fio = 
100 meV. The insets indicated by an arrow, show the structure 
of the solid blue curve in comparison with the input electron- 
phonon spectral function (red dashed curve). See text for 
further discussion. The lower frame is also for the Lorentzian 
distribution used for the upper frame but now = 500 meV. 



size and shape for each doping value. 

As we noted with reference to the top frame of Fig. [U 
the logarithmic singularities at tu = ±uje in Eq. ([6]) for 
Rel](ci;) no longer appear when the chemical potential 
/iQ = and the bands are filled exactly to the Dirac point. 
This is shown in Fig. [S] (top frame). Only a small "kink" 
(i.e., change in slope) remains at lu — ±lue- This fact 
has important consequences for experiment. While the 
density of states as a function of lu still knows about the 
renormalization factor of 1 + (middle frame) there 
is no additional logarithmic-type signature of the Ein- 
stein phonon mode. Consequently, compared with the 
examples of Fig. [51 only very small peaks are seen at 
LU = ±lue in the derivative, as shown in Fig. [5] (bottom 



11 



200 I I I I I I I I I I I I I I I I I I I I 




-200 I I I I I I I I I I I I I I I I 




cj(eV) 

FIG. 9: (Color online) Top frame: Self-energy T.{uj) for no — 
in units of meV for the real (solid curve) and imaginary part 
(dashed curve) . The density of states for this case is shown in 
the middle frame in units of meV with the bare DOS given as 
the dotted curve. Derivative of the DOS normalized to A'o is 
shown in the bottom frame with the dashed black line is the 
bare case. In all cases, the solid blue curves includes phonons. 



frame). Nevertheless, one can see the 1 + A®*^ renormal- 
ization at low frequency in this frame, given as the mag- 
nitude of the height of the curve near zero frequency. In 
their scanning tunneling microscope (STM) results, Li et 
al;^ find structure in the DOS at 155 meV. In addition, 
they find a renormalization which they determine to be 
A''*^ = 0.26, for a case of doping very close to the charge 
neutrality point. This case was addressed in our previ- 
ous short communication^^ which was aimed towards a 
detailed evaluation of the data analysis procedure used 
by Li et al. for extracting a value of the mass enhance- 
ment. From that we conclude that the data may support 
a larger value of A, possibly as large as 0.4. 



IV. ANGLE-RESOLVED PHOTOEMISSION 
(ARPES) 
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FIG. 10: (Color online) Renormalized energies (solid blue 
curves) i?*, as a function of k in units of eV/{hvo) for a case 
with fio — 400 meV. Twice the bare and dressed chemical 
potential are indicated by vertical arrows, long and short, re- 
spectively. The bare curves are shown as dashed black lines 
and the dotted is a line chosen to fit the dressed curve asymp- 
totically at large negative energy. 

The charge carrier spectral density A{\i,uj) given by 
Eq. ID) can be measured in angle-resolved photoemission 
experiments (ARPES). The technique is directional and 
measures A(k, lu) as a function of lo for any direction and 
magnitude of momentum k. Here, because of the symme- 
try only e = ±?iuo|k| enters. As we have already noted, if 
the imaginary part of the electron-phonon self-energy is 
infinitesimal, A{zLe, ui) would reduce to a delta function 
and provide the dressed energy Ek corresponding to the 
bare through the equation Ek — KeT,{Ek) + fi — Ck = Q- 
In Fig. [TOl we show results for Ek in the case of /io = 400 
meV. The dashed black lines are the bare dispersions and 
the solid blue curves, the renormalized energies. Several 
features of this figure need to be noted. The renormal- 
ized dispersions show a clear "kink" at the peak phonon 
energy ujq — 200 meV. In the calculations, we used our 
Lorentzian model for the distribution of phonon energies 
about luq otherwise the "kink" would actually show a 
logarithmic singularity. For a simple metal with a con- 
stant DOS around the Fermi energy, the same phonon 
anomaly would arise, but at higher energies (below the 
Fermi energy) the renormalized curve would rapidly re- 
turn to the bare value as the real part of 5](w) drops to 
zero. This is not true for graphene because the density 
of electron states increases linearly with increasing en- 
ergy and this means that |ReI](cij)| increases as we saw 
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in the top frame of Fig. [T] The phonon structures are 
superimposed on this more gradual increase in |ReI](ct;)| 
versus w. The dotted red Une in Fig. [10] brings out this 
feature clearly. It has been drawn to match the dressed 
dispersions at large negative energies. In contrast to what 
would be the case for the bare dispersions, this asymp- 
totic line when extended towards the Fermi energy does 
not cross through the Fermi energy. This has been noted 
in the ARPES data of Bostwick et al.^^. In our calcu- 
lations this is a direct consequence of the fact that the 
electron-phonon renormalizations continue to increase at 
high energies because of the continuing increase of the 
bare density of states. In the actual experimental data, 
other renormalizations would also contribute. While the 
slope of the bare dispersion curves (dashed black lines) 
gives directly the bare velocity vq the slope of the asymp- 
totic dotted line is smaller and would give a smaller value 
of the effective velocity of light for graphene should it 
be used to define the bare bands. It is close but not 
the same as the renormalized effective velocity of light 
t;o/(l + A''*^). The vertical arrows in Fig. [TU] identify 
renormalized (shorter blue arrow) and bare (longer black 
arrow) chemical potentials, respectively. It is clear that 
/X is about 15% smaller. This is close but does not quite 
correspond to the factor of 1 + — 1.19 in our calcula- 
tions indicated in Eq. ([^0]) which is strictly valid only for 
/io <C Wo- Nevertheless, the differences are small. Here 
the factor of 1.19 gives exactly the slope of out of the 
Fermi energy. 




-500 -400 -300 -200 -100 □ ICQ 200 300 400 500 



k 

I 



FIG. 11: (Color online) Color map of interacting dispersions 
Ek in units of meV as a function of k in units of meV/{hvo). 
The bare chemical potential /uo = 400 meV and a Lorentzian 
phonon spectrum with ljo = 200 meV were used. 



Of course, the imaginary part of S(ti-') in Eq. ([T]) for 
Lu ^ would never be zero and we would not have per- 
fectly well-defined dressed dispersion curves. In such a 
circumstance, it is useful to use a color map to show 
the constant energy contours of A(e, w). This is done in 
Fig. [TT] Cuts at different values correspond to different 
colors as indicated on the scale shown on the right (units 
of meV~^). Because we have used a phonon spectrum 
which starts at = 170 meV for this figure {ujq — Sc), 
there is no quasiparticle lifetime above —170 meV. Be- 
low this energy the curves widen and reflect the broad- 
ening. We see clearly the peak of the phonon anomaly 
at —200 meV but the identification of the position of the 
Dirac point becomes somewhat ambiguous because of the 
broadening of the curves. The renormalized \fj,\ — 344 
meV for this case and according to Fig. [7] this almost co- 
incides with \ujd\- The pinching at —600 meV is due to 
the small value of the imaginary part of the self-energy at 
—fiQ—oJo as we have already seen in the middle frame of 
Fig. [T] for a different value of chemical potential. Other 
interactions such as electron-hole particle pair formation 
or coupling due to plasmons are expected to lift this 
pinching and provide a finite lifetime. It is clear from 
this figure that the constructions used in Fig. [10] to iden- 
tify the dotted asymptotic line as well as the interacting 
chemical potential are not as well-defined when broaden- 
ing is included in the theory. 



V. ELECTRON-PHONON INTERACTION ON 
DC CONDUCTIVITY 

Next we turn to the optical conductivity. In units 
of the universal background conductivity o-q = 7re^/2/i, 
where e is the electron charge and h is Planck's con- 
stant, the real part of the frequency-dependent conduc- 
tivity a{T, f2) at frequency and temperature T is given 
by 



^ ' = K / duj[f{u)- f{uj + VL)] ede[A{e,uj)+A{~e,uo)][A{e,uo + n) + A{-e,uj + n)], (31) 

CTO " J-oo Jo 

I 

where we have ignored vertex corrections. In metal extra weighting of (1 — cos9) in the calculation of the 
physics, this can be incorporated approximately by an scattering rate changing it from a quasiparticle to trans- 
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port rateJJi This has recently been verified to hold as 
well for the specific case of the DC conductivity of 
graphene by Cappelluti and Benfatto^. The function 
f{u!) = l/[exp{f]uj) + 1] is the Fermi occupation factor 
with /?, the inverse temperature. We start with a general 
observation that 



(32) 



To establish this symmetry, we note first that for negative 
values of the chemical potential /i, the carrier spectral 
function satisfies 



A 



where we have used Eqs. ([7]) and 
be rewritten as 



(33) 

Also, Eq. (I3l|) can 



a^<o(T,r!) ^ 4 ^^[j(^) „ f^)] ede[A{e,-cj) + A{^e,-u;)][A{e,-iJ^n) + A{-e,-u;-n)]. (34) 

C^O ''U-oc Jo 
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Noting Eq. then Eq. ^ follows. 

Taking = in Eq. ([3T|) gives the DC conductivity at 
finite T, 



J 



o-Dc{T) 



-if du;^^ f " ede[A^{e,uj)+A^{-e,uj) + 2A{e,uj)A{-e,uj)]. 



OO 



(35) 



r 



We have verified numerically that in Eq. (|35p . we can simplifies the result and we get 
send Wc to infinity within good approximation. This 



J 



(^Dc{T) 



duj 



dm 

duj 



r 



tan 



a(ti;) 

rv) 



(36) 



with r(a;) = — Iml](a;) and a(cj) = cj — ReI](cL;) + /i. Note 
that for S(w), we have neglected the temperature depen- 
dence that would enter Eq. Q at finite temperature'^^ 



and hence our results here are for low T. At zero tem- 
perature —df{u!)/duj in Eq. ([55)1 becomes a Dirac delta 
function centered around lo — (i.e., S(lu)) and 



J 



<TDc{T = 0) 



fi - ReS](u; = 0) 
T{uj = 0) 



r(cj = 0) 



- Rel](w = 0) 



tan 



_^ f fi - ReT.{uj = 0) 

r(tj = 0) 



(37) 



r 



But /i — Rel](w = 0) in ([37| can be replaced by /iq, the 
bare band chemical potential from Eq. pT|) . Also, at 
= and T — 0, the electron-phonon interaction does 
not contribute to the scattering rate which therefore re- 
duces to the residual scattering rate rj, i.e., r(a; = 0) = ry, 
which is small but finite for graphene. Thus, the electron- 



phonon interaction has entirely dropped out of Eq. (|37p . 
which reduces to its bare band value. This means that 
the familiar resul^iii that the electron-phonon interac- 
tion does not change the DC conductivity of ordinary 
metals also holds for graphene. 
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For fi - ReE(w = 0) = /io < in Eq. ((371), we get the 
universal DC limit )^^''^^ unrenormalized by the electron- 
phonon interaction 



8 4e^ 

CTDciT = 0) = — CTo = ^ 



(38) 



which differs from the universal AC background value 
(To by a factor of S/ir^ ~ 0.81. The first correction for 
finite /iQ is included when Eq. psp is multiplied by (1 + 
(/io/??)^/2). In the opposite limit of /i — RcE(cj = 0) ^ 77, 
we get instead 



<TDc{T = 0) = 



- Rel](w 0)] 
2^ 



h 2r] ' 



(39) 



Here, 2/io is the optical spectral weight removed from the 
universal background at finite chemical potential (in the 
bare band) which reappears as a Drude type contribu- 
tion about = 0. Note that l/r = 277 is the transport 
scattering rate. In this limit, the DC conductivity is 
no longer universal, i.e., rj does not drop out as it did 
in Eq. psp , however the electron-phonon interaction has 
dropped out. This is similar to the well-known resuHjiSiiii 
in conventional metals, that the electron-phonon interac- 
tion does not change the value of the DC conductivity. 
This remains true for graphene at T — 0, and can be 
traced to Eq. (|57|) . <7dc{T = 0) does not depend on the 
electron-phonon interaction because of the appearance 
of the chemical potential shift factor ReS(u; = 0), which 
changes /j, to /io, the bare band value, and r(aj = 0) 
reduces to the residual scattering rate. 

Returning to Eq. (j35p . it is useful to separate intra- 
band and interband contributions. The first two terms 
are intraband because they involve transitions within 
the same Dirac cone while the last term involves tran- 
sitions between lower and upper cones. In a conven- 
tional metal only the intraband contribution arises in 
discussions of the DC conductivity. In Fig. [121 we 
show results for uociT — 0) in units of ctq as a func- 
tion of (^ — Rel](a; = O))/?? = ^Jio/v■ The sofid (red) 
curve is the interband piece, the long-dashed (blue) curve 
the intraband and the short-dashed (black), the sum. 
Both the intraband and interband transitions contribute 
equally to the universal DC value which corresponds to 
[ji — Rel](a' — 0))/ry = /^o/?? = 0, i.e., the lower cone 
is filled to the Dirac point and the upper cone is empty. 
As /io is increased and the material is doped away from 
the neutrality point, the interband contribution decreases 
while at the same time the sum increases rather substan- 
tially. What determines the scale of this increase is the 
value of the ratio of the chemical potential to the resid- 
ual scattering rate rj. While the value of 77 is somewhat 



uncertain for graphene and will vary with sample quality, 
it is expected to be small and of order of 1.0 meV so that 
the value of the DC conductivity can be increased a lot 
for easily attained values of the chemical potential, say 
/io of order a few hundred mcV. 

Returning to Ea. ((36| which gives the DC conductiv- 
ity at finite temperature T, we note that the value of 
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FIG. 12: (Color online) DC conductivity a(Sl = 0) at T = 
in units of the universal background value ao — -Ke^ /2h. The 
solid red curve gives the interband contribution, the long- 
dashed blue curve, the intraband and the short-dashed black 
curve, the sum as a function of [/i — ReE(aj = O)]/?; = 
/io/r/, with /i the chemical potential, /to its bare band value, 
ReE(ti; = 0) its shift, and r/ the residual quasiparticle scatter- 
ing rate. 



Rel](w) for finite lo now enters the formula even for /i = 
(Dirac point is at the Fermi surface) and consequently 
the DC conductivity is affected by the electron-phonon 
interaction. It is only for T = 0, /fo = that wc get 
ubc = 0.81(To independent of electron-phonon renormal- 
ization. 

In a conventional metal, if only impurity (elastic) scat- 
tering is accounted for, the DC conductivity does not 
change with increasing temperature. To get a change it 
is necessary to include inelastic processes. This is not so 
for graphene. Including only residual scattering and no 
electron-phonon renormalization, the formula for the DC 
conductivity Eq. reduces to 
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In Fig. [131 we show on a logarithmic scale a{^l = 0, T)/aQ 
as a function of /i = fi/T for various values of 77 = rj/T. 
Here, = /iq, the bare band chemical potential, as T,{uj) 
has been taken to be zero. We first note that a,t p, — 0, 
the universal limit is no longer universal unless the tem- 
perature is much smaller than the residual quasiparti- 
cle scattering rate 77. An analytic expression valid for 
T/r] < 1 is^ 



4e2 
irn 



9 u 



(41) 



For rj — T, the DC conductivity has increased by more 
than a factor of two over its T = value. It is clear 
that in clean systems, it is necessary to go to very low 
temperatures in order to observe the universal limit. It 
is also clear that the chemical potential is to be small 
compared with rj. It is the residual scattering rate which 
sets the scale on T and /i. The conditions to observe the 
universal limit are T <^ rj and fj, <^ r]. We turn next to 
finite O in the microwave or Terahertz region where is 
of order 77, and study the evolution of the conductivity 
from its DC value to its universal AC background value. 
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FIG. 13: (Color online) The DC conductivity g{VL = 0, T) at 
finite, but small, temperature T in units of ao — ne^ /2h as a 
function of i-io/T for several values of ry/T with fio, the bare 
chemical potential and rj, the residual quasiparticle scattering 
rate. Here, for simplicity, phonon renormalizations are taken 
to be zero. 



VI. LOW FREQUENCY CONDUCTIVITY 



In Fig. [Ml we show results for a{n) at T = in the 
frequency region spanning a few times the energy of the 
quasiparticle scattering rate 77. The top left hand frame 
shows a{fl)/ao versus 51 up to 100 meV for zero chem- 
ical potential with rj — 1 meV while the lower frame is 
the same but now rj — 2.5 meV. In both cases the long 
dashed (red) curve is the interband contribution and the 
short-dashed (blue) is the intraband. The sum is the solid 
black curve. We see that at = 0, the universal limit 
is reached for both values of r] as expected. As ft is in- 
creased a{Q) increases toward the universal background 
value (To = 7re^/(2/i). The energy scale on which satu- 
ration is reached is set by the value of rj and hence will 
vary with sample quality. The two right hand panels are 
similar but are for finite chemical potential, with fiQ — rj 
for the top and fiQ = brj for the bottom, and rj = 2.5 
meV in both cases. It is clear that the conductivity in 
the low region is greatly affected by finite /xg and is 
non-universal. This arises because at finite yuo, the in- 
terband transitions below 2fio are Pauli blocked and the 
optical spectral weight involved is transferred to the in- 
traband transitions. This provides a Drude peak with 
finite effective plasma frequency proportional to 2/io. In 
the top right panel of Fig. [141 A*o = so that we are not 
in the universal DC limit regime and its value is in fact 
already larger than the universal background value (Tq. 
As n is increased, a(ri)/ao drops, shows a minimum be- 
fore rising again to reach its universal background value. 
The energy scale for this final rise to saturation remains 
ij as in the left hand panels but this is because 77 is of 
order /xg. For a case iiq > i] as in the lower right panel, 
specifically with /io = IO77, the Drude at small f2 is much 
more pronounced as compared with the upper frame and 
it is now 2/i which sets the energy scale for recovery of 
the conductivity to its universal background value. An 
interesting feature of the interband contribution to the 
conductivity (long-dashed red curve) in this case is the 
long tail extending to fl = which is nearly constant 
before its main rise around Q = 2/x. This behavior can 
be understood analytically as we now show. 



At zero temperature, for small but finite Q the con- 
ductivity is given approximately by 



J-n "^^Jo ^ 



ede 
2" 



[cj(l + Aoff)±e + Aio]2+77^ 



[{uj + + Acff) ± e + Aio]2 + 



(42) 



where we have noted the lo is also small and have used 



— \cSUJ for ReS(ci;), after subtraction of the constant 
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piece. The imaginary part of the self-energy due to 
phonons is zero so that only the residual contribution re- 
mains. Should the acoustic phonons have been included 
in the model, ImS](a;) would be finite, but this effect is 
small in graphenc and neglected here. The sum over ± 
leads to four terms in Eq. (f42|) . two intraband contri- 
butions (same sign) and two interband (opposite signs 
of e in pair of terms). The integral over energy can be 
performed analytically and the sum over -t- and — , i.e., 
upper and lower Dirac cone performed. Starting first 
with the intraband piece and writing f2 = 0(1 -I- A°^) 
and U! = uj{l + X"^), we get 



FIG. 14: (Color online) The conductivity a{Q.) at T = in 
units of ao = Tie^ /2h as a function of Q emphasizing the low 
frequency region. Left frames are for zero chemical potential 
{fio = 0) and quasiparticle residual scattering rate rj — 1.0 
and 2.5 meV. Right frames are for rj = 2.5 meV with — rj 
and lOr;. Long-dashed red curve is the interband contribu- 
tion, short-dashed blue is the intraband and solid black is the 
sum. Electron-phonon renormalizations are included using 
the parameters previously discussed. 
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(43) 



Assuming u!,Q,,ri ^ /io, expression ([43]) greatly reduces 
and we get 

0'»«tra(^) ^ Mf^ol 

cro TT 02(1 + A<=ff)2+ 4772 ^ 

for the intraband contribution to the conductivity when 
O <C /io; as well as, ^ /ig. The form is a Drude 
with effective optical scattering rate 277/(1 + X'^^) and 
effective plasma frequency of 4/io/[7r(l -t- A"**")]. While 
we have made approximations to derive Eq. (|44p . we 
have verified in numerical work for Fig. [14] right bottom 
frame that the long-dashed blue curve follows ([33]) very 
well. The result ([33]) is the same as found in conven- 
tional metals. The coherent part of the conductivity is 
a Drude form with scattering rate and plasma frequency 
renormalized by (1 -t- A''^) and the DC limit remains un- 
renormalized by A°^. Here for simplicity, we have used 



a constant residual scattering rate. For graphene, this 
scattering rate can itself be frequency dependent and this 
can change the Drude line shape which then reflects this 
additional energy dependence as described in Rcf.^-. 

Application of the same sort of algebra for the inter- 
band contribution to the conductivity in the limit O and 
77 /xq, gives a leading contribution which is frequency 
independent: 

(Jinteri^) 1 ^fiQ-T] 

= 5 5", (45) 

(To TT 772 + 

(O < uje) and in sharp contrast to Eq. ([33]) is indepen- 
dent of the electron-phonon renormalization parameter 
X°^ . Thus, while the intraband Drude contribution is 
renormalized by the electron-phonon interaction in the 
same way as in ordinary metals, the interband is not. 
Equally remarkable is that ([¥5]) . in contrast to ([H]) . is 
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FIG. 15; (Color online) The interband conductivity (Tenter (fi) 
at T = in units of ao — 7re^/2/i as a function of photon 
energy Q emphasizing the low frequency region, for varying 
values of Ho/Vi with rj — 2.5 meV. Electron-phonon renormal- 
izations are included as described previously. 
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FIG. 16: (Color online) AC conductivity a{Q,)/ao for two 
values of chemical potential /io = 150 meV (top frame) and 
fiQ = 400 meV (bottom frame), one below and one above 
the phonon energy of cue = 200 meV. The dashed red curve 
Wcis obtained without electron-phonon renormalization and 
the solid blue is with. 



VII. INFRARED CONDUCTIVITY 



independent of photon energy, as we already noted in the 
lower right hand frame of Fig. [TJ] where we present nu- 
merical results. The interband contribution (long-dashed 
red curve) remains essentially constant at its DC value 
even for ^ 2/i at which point it begins a steep rise and 
eventually makes the major contribution to the universal 
background value. At small 51 it is the reverse that holds, 
the interband piece is small while the intraband Drude- 
like contribution is dominant. Additional results for the 
interband contribution to the conductivity for photon en- 
ergies less than and slightly above the value of twice the 
chemical potential are shown in Fig. 1151 Here we span 
values of from to 10 times 77, as labeled on the figure. 
Only the first curve (solid black) satisfies the condition 
for the universal DC limit and in this case, the curve 
starts at a value of 4/7r^ ^ 0.4 and makes half the contri- 
bution to the universal DC conductivity. Furthermore, 
it rapidly rises to its universal AC background value on 
the scale of 27], with rj = 2.5 meV. All other curves start 
with (T{n)/(To below 4/7r^ at = 0. This is true even 
for the red dotted curve for which /iq = ??• The last 
few conform well to the analytic prediction of Eq. ([l5|) 
that ainter{^)/<yo = (2a<o/7I'?7)/[1 + ilJ'o/ri)'^], which we 
see extends to ^ 2/i and, as noted before, the main 
rise towards the universal background value is set by 2/i 
rather than by rj. The few intermediate curves show the 
evolution from one regime to the other. 



So far we have examined the small photon energy range 
of the AC conductivity as compared with the value of the 
phonon energy lue = 200 meV in our model. In Fig. 1161 
we show results to $7 = 1 eV for fio = 150 meV (top 
frame) and ^0 = 400 meV (bottom frame). The dashed 
red curve, which is included for comparison, give results 
without the electron-phonon interaction (bare band case) 
but with a phenomcnological impurity transport scatter- 
ing rate of 2rj with rj = 2.5 mcV (constant). The solid 
blue curve has phonons. The first thing to note when 
comparing these curves is that the main rise in a(fl)/ao, 
indicating the increase in interband transitions, occurs 
at twice the value of the dressed chemical potential in 
the solid curve while it occurs at 2/io for the bare bands. 
This translates into a considerable shift downward of this 
prominent threshold for absorption. This is easily under- 
stood with the help of Fig. [TUlwhere bare band and renor- 
malized dispersion curves are shown. Absorption of light 
for the interband case proceeds through vertical transi- 
tions from an occupied state in the lower Dirac cone to 
an unoccupied state in the upper Dirac cone, since the 
photon transfers no momentum to the electronic system. 
This is shown as the arrows which applies to the lowest 
energy transition allowed by the Pauli exclusion princi- 
ple, which blocks transitions below this energy. For the 
bare clean bands this energy is just 2/io (large black ar- 
row), twice the bare chemical potential. But in the in- 
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teracting band the first transition occurs rather at 2/i 
as shown on the left hand side of the figure by a blue 
vertical arrow which is clearly shorter in length than is 
the black arrow. In reality, of course, there are finite 
lifetime effects which broaden the initial and final states 
(except for the one right at the Fermi energy) and this 
is reflected in a rounding of the absorption edge. In our 
numerical work for bare bands, we have included a small 
residual scattering rate on our states so that even for the 
red dashed curve a small rounding of the interband edge 
is seen in the figure. Note that this shift is also present 
in the previous Figs. [HI and [T5l 
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FIG. 17: (Color online) AC conductivity a{Q)/ao showing 
separately intraband (short-dashed blue curve) and interband 
(long-dashed red curve). The solid black curve is the sum of 
the two. The chemical potential /io — 400 meV. 

Another feature to be noted when considering Fig. 1161 
is the contrast between top and bottom frame for which 
fiQ goes from a value of 150 meV, which is less than 
uje = 200 meV, to a value of /zq = 400 meV which is 
larger. The large value of /xq allows the Holstein phonon- 
assisted sideband to be prominently revealed. It starts 
at — 200 meV in the blue curve, while for the bare 
band it is not present. What is important to realize is 
that both interband and intraband processes contribute 
to this absorption. The intraband contribution is easily 
understood from conventional metal theory and arises 
when a photon not only induces an electronic transition, 
i.e. creates an excited hole-particle pair, but also creates 
a phonon in the final state. But there is also an interband 
contribution to the Holstein sideband as can be seen in 
Fig-im where the total conductivity (sohd black curve) is 
decomposed into an intraband (short-dashed blue curve) 
and an interband (long-dashed red curve) contribution. 
To understand the behavior of the interband curve we 
refer to Fig. I18[ where we show results for the probabil- 



ity of occupation of a state k, nk, which can easily be 
calculated from the spectral density ^(k, cj) of Eq. 
At finite temperature 

/oo 
f{uj)Aik,uj)dLj. (46) 
-oo 

If, for the moment, we ignore damping effects, the elec- 
tron spectral density ^(k, lj) becomes a delta function 
and the probability of occupation of a state (e^ — /io) 
in the bare band becomes nk = f(hvo{k — kp)) for the 
Ek = +hvok branch and nk = f{—hvo{k + kp)) for the 
Ek = —hvok branch, which is just the thermal occupation 
function. When the EPI is included, an interesting ana- 
lytical result can still be obtained if we limit ourselves to 
the region of small energies away from the Fermi surface. 
In that case, we can approximate Rel](a;) by Eq. (jlSp 
and obtain 

"k^TTA^^^^^'o^^"^^^^' ^^^^ 

where Vq (defined below Eq. HI])) is the renormalized 
effective velocity of light, kp is the bare band Fermi mo- 
mentum and the equation is valid only for (k — kp) small. 
We can immediately see from (|47p that the jump in oc- 
cupation at the Fermi surface in the interacting system 
at r = has been reduced below one by the factor 
1/(1 + A''^). Our numerical results with and without 
phonons are shown in Fig. 1181 The dashed blue curve 
is close to unity except very near the Fermi surface at 
k = kp in our notation. We have used a residual scatter- 
ing 7] = 0.001 meV. When phonons are included, we ob- 
tain the solid red curve which also shows a sharp drop at 
the Fermi surface, but the jump is now reduced below one 
by a factor, 1/(1 + as we saw in Eq. (|T7)) . In addi- 
tion, this curve has very significant tails beyond the Fermi 
level and at negative energies below this point, the prob- 
ability of occupation remains reduced from the value of 1 
on the scale of eV reflecting the fact that phonon renor- 
malizations persist to high energies in graphene as we 
saw in Fig.[TJ Thus since nk is always less than 1, Pauli 
blocking is lifted on interband transitions. The initial 
state in the occupied lower Dirac cone of Fig. [TU] remains 
occupied with finite probability but at the same time 
a final state below the Fermi level can now accommo- 
date a photo-excited electron since this state is not occu- 
pied with probability one. This is illustrated by the inset 
schematic of Fig. [181 where the large black arrow repre- 
sents the case where the particle is promoted to above the 
Fermi level, which would be a transition with minimum 
energy 2^. However, the finite probability for holes to 
exist below the Fermi level gives rise to transitions such 
as shown by the short red arrow, and interband absorp- 
tion can now proceed for energies which are less that 2fi. 
This partial lifting of Pauli blocking is responsible for the 
long tails extending to O = seen in the lower right hand 
frame of Fig. [HI which allow the interband transitions to 
contribute to the DC conductivity. They also allow for 
additional phonon-assisted interband transitions to start 



19 



at r2 > as seen in the long-dashed red curve of Fig.lTTl 
An interesting point to be aware of is that for large values 

1.2 I I , , , , J 




0.0 ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' 
-1 -0.5 0.5 

Ek-Mo (eV) 



FIG. 18: (Color online) Probability of occupation of state 
Efc for two cases: bare band (dashed blue curve) and with 
electron-phonon interaction included (solid red curve) . A very 
small residual scattering of rj — 0.001 eV is also present in 
both cases. The inset is a schematic which illustrates the 
renormalized energy bands filled to the Fermi level Ef with 
finite probability for some holes to exist below the Fermi level. 
Interband transitions are now possible for energies below 2/i. 

of the chemical potential, a regime in which a Drude and 
its boson-assisted side bands is clearly revealed, the side- 
bands will not have the same relationship to the Drude 
as they would in a simple metal. We expect an additional 
contribution in this energy range from interband transi- 
tions, as shown in Fig. [171 We saw in Eq. (|47p that the 
jump at k — kp in the occupation probability is renor- 
malized by the EPI by a factor of 1/(1 -I- X°^). We also 
saw in Fig. [18] that these renormalizations persist to high 
energy away from the Fermi surface and the occupation 
factor Uk is significantly reduced below one allowing for 
the relaxation of Pauli blocking on interband transitions. 
On the other hand, these renormalizations drop out of 
the value of the interband conductivity in Eq. (j^J) valid 
for small photon energy and Fig. [T5l shows this to be true 
for a considerable range of energy above = 0. This is 
a clear illustration that the EPI can renormalize various 
quantities in quite different ways. 

Another interesting aspect of the redistribution of opti- 
cal spectral weight that is brought about by the electron- 
phonon interaction is illustrated in Fig. [THlwhere we show 
the results for the partial optical sum denoted by 
and defined as 

7(0) = ^^dn', (48) 
with variable upper limit Q on the integral. The red 



solid curve is the bare band case but with a small resid- 
ual quasiparticle scattering rate rj — 2.5 meV included 
in the numerical evaluation of Eq. (|3ip . In this case, 
the entire spectral weight 2/io lost through Pauli block- 
ing of the interband transition is found in the Drude. 

rises sharply to a height of 0.8 eV on an energy 
scale given by 2?], then remains constant until = 0.8 
eV = 2/xo where it starts rising again in a linear fashion. 
The solid black curve which includes phonons is more in- 
teresting and is to be compared with the red curve. It 
too rises out of zero very rapidly on the scale of 77 but 
does not go all the way up to a height of 2/Lto rather it 
shows a saturation plateau at 2/xo/(l -I- X^^) = 2/i in the 
region below ue — 200 meV. At uje, intra and interband 
Holstein processes start to kick in and begins to 

rise slowly until ^ 2/i is reached where a considerable 
change in slope occurs because the major contribution of 
interband transitions starts coming in (as can be seen by 
comparing the blue short-dashed intraband curve with 
the solid black curve for the total amount). By oj — 2/io, 
the black and red curves have effectively merged and 
the optical spectral weight redistribution brought about 
by the EPI is in balance at this point. At larger ener- 
gies fi ~ 2 eV, the solid black curve which is based on 
the renormalized conductivity remains below the solid 
red curve. This is because the electron-phonon interac- 
tion has reduced the value of the universal background 
slightly below ctq — ne^ /2h (also seen in Fig. [16] at large 
O, where the solid blue curve is slightly below the red 
dashed curve). This is understandable. We saw in Fig. [S] 
that the electron-phonon interaction has a profound ef- 
fect on the band structure in the energy region around 
the band edge. In this region, the DOS is considerably 
depleted below its noninteracting value and to conserve 
states tails appear beyond the bare cut off Wc- Thus 
in optical experiments, spectral weight is removed below 
the bare optical cut off which is transferred to higher en- 
ergies. The slight reduction below one of the universal 
background is generic and we would expect it to be a 
feature of interactions in general. 



VIII. SUMMARY AND CONCLUSIONS 

It is well known^°'^^ that in wide band metals for 
which the band density of states (DOS) is essentially 
constant on the scale of a few times a phonon energy, 
the electron-phonon interaction drops out of the renor- 
malized DOS. Because of the Dirac nature of electron 
dynamics in graphene the bare band DOS is linear in en- 
ergy rather than constant and consequently an image of 
the phonons is retained in its dressed DOS and a first 
derivative of N{lo) provides an ideal baseline to study 
boson structures. For a bare Dirac band this derivative 
gives minus or plus one with the change in sign occurring 
at — /ioj i-S- at minus the bare chemical potential which 
also coincides with the Dirac point. Several modifications 
arise when the electron-phonon interaction is included. 
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FIG. 19: (Color online) The partial optical sum /(fi) defined 
in Eq. (|48p vs S7 (variable upper limit on the integral). The 
red long-dashed curve is for the bare band case with r\ = 2.5 
meV and is for comparison. The solid black curve has phonon 
renormalizations. The short-dashed blue curve includes only 
the intraband contribution in the presence of phonon renor- 
malizations. 



First, the Dirac point position is shifted and this leads 
to a change in the energy of the jump from minus to 
plus one in the first derivative. This jump remains sharp 
provided the chemical potential is in magnitude smaller 
than the phonon energy. But when the electron-phonon 
lifetime is finite at the Dirac point, the jump becomes 
smeared. In addition there are structures at icj^; for 
an Einstein spectrum. When a distributed phonon spec- 
trum is used, the structure in dN(uj)ldLo reflects not only 
its magnitude but also it shape. With increasing doping 
or more precisely magnitude of chemical potential, the 
phonon structures are enhanced because the Fermi level 
falls in a region of higher density of states and conse- 
quently of larger strength of the effective electron-phonon 
interaction. A point of note is that at the Dirac point, 
the explicit logarithmic phonon structures of the finite 
doping case disappear, although the DOS still retains 
its (1-1- A"*^) mass renormalization. Small features re- 
main in the derivative and there is some evidence from 
STM2^ for both of these effects in this limit. The method 
should be applicable to other scattering mechanisms such 
as electron-hole and plasma excitations. 

Another relevant observation is that in the DOS, we 
find that for /ip > w_e, the interactions lift the Dirac point 
from a value of zero in the bare case to a finite value, and 
that the DOS rises quadratically rather than linearly out 
of this point. This is due to a finite scattering rate at 
the Dirac point. Recent STM experiments'*^ appear to 
see this feature of lifting and accompanying quadratic 
behavior. Further investigation of this point would be 



warranted. 

We have considered separately the effect of the 
electron-phonon interaction (EPI) on the intra and inter- 
band contributions to the optical conductivity. We find 
that the DC conductivity remains unrenormalized with 
equal contributions from both processes when the chem- 
ical potential is zero. With increasing doping through 
application of a gate voltage in a field effect device or 
by seeding the graphene surface with potassium or other 
atoms, the interband contribution decreases and becomes 
negligible when the transport scattering rate r] becomes 
much smaller than the chemical potential. The con- 
ditions for the observation of an universal DC limit is 
fj. <^ rj, as well as, T <^ rj. For T > -q, there is a rapid 
increase in DC conductivity in contrast to ordinary band 
metals for which it is independent of temperature. For 
zero doping, the scale on which the universal DC value 
Ae'^/nh rises with increasing photon energy to its univer- 
sal background value of 7re^/2/i is set by the transport 
scattering rate r/. By contrast, but in complete paral- 
lel to what is known to apply in simple band metals, 
as the chemical potential is increased, a clear AC intra- 
band Drude conductivity is revealed which is renormal- 
ized in two ways. Through analytic techniques, we show 
that the transport scattering rate l/r is to be replaced 
by 1/[t(1 -I- A°*)] and the plasma frequency squared re- 
duced to + X"^). Here, A^^ is the electron ef- 
fective mass renormalization at the Fermi energy which 
varies with doping. In addition, optical spectral weight 
is transferred to higher energy in the form of a Holstein 
phonon-assisted absorption sideband. Also, interband 
transitions, not part of the theory of simple metals, pro- 
vide additional absorption. At small frequencies below 
twice the value of the chemical potential (^) there is a 
small but finite nearly frequency independent interband 
absorption which shows an additional Holstein side band 
for the case 2fj, > uje (the phonon energy), with an ad- 
ditional rapid increase towards its universal background 
value tTo = 7re^/2/i at 2/i. The near constant in- 
terband absorption below 2/i can be traced to a lifting 
of Pauli blocking brought about by the EPI. The inter- 
actions decrease the probability of occupation of a plane 
wave state below the Fermi energy to a value less than one 
so that such states can still be used as final states for in- 
terband transitions although these are greatly reduced in 
optical spectral weight. In pure graphene, the increase in 
conductivity at f2 = 2/io represents a vertical jump from 
to (To- When the EPI is included not only does the 
chemical potential shift to a smaller value but the jump 
at 2^ is also smeared due to finite lifetime effects. While 
at small frequencies in the Drude region the conductivity 
can exceed its background value ctq, for energies greater 
than 2/i it is always slightly smaller than ctq- Consider- 
ation of the partial optical sum rule shows that missing 
optical spectral weight at high energy is transferred to 
even higher energies as compared to the bare band case 
because the EPI provides states beyond the bare band 
cutoff on the scale of the phonon energy. It also shows 
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that at small energies the readjustment of optical spec- 
tral weight between bare and interacting case is in near 
balance at Q ~ 2/io- 
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